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Using a recently developed approach to constructing covalent network structures from linear carbyne, we 
identify by ab initio calculations a new carbon allotrope in R3m (Dl d ) symmetry that comprises six-fold 
helical chains with alternating spt-type single and double bonds along the chains that are connected via 
zigzag benzene rings. This 6-fold carbene is characterized as a three-dimensional three-connected chiral 
crystalline modification of graphite. Phonon and electronic band calculations indicate that this new 
structure is dynamically stable and is a semiconductor with a band gap of 0.47 eV, in contrast to the 
semimetallic nature of graphite. Simulated x-ray diffraction patterns of the 6-fold carbene provide an 
excellent match to the previously unexplained distinct diffraction peak of a new carbon allotrope found in 
recent detonation experiments. These results establish a new carbon phase and offer insights into its 
outstanding structural and electronic properties. 



Carbon's extremely versatile bonding ability produces a rich variety of structural allotropes, among which 
graphite and diamond have long been known for their outstanding properties. All carbon atoms in 
graphite are connected via the sp 2 hybridization bonding with a bond length of 0.142 nm and bond angle 
of 120°, and the graphite honeycomb lattice can be viewed as a planar molecule comprising benzene rings or 2- 
fold polyacetylene-like zigzag chains 1 . The strong in-plane aromatic n- conjugation makes the two-dimensional 
three-connected (2D3C) graphite the most stable allotropic form of carbon. In a diamond single crystal, all carbon 
atoms are connected via the sp 3 hybridization bonding with the tetrahedral bond angle of 109.47° and a bond 
length of 0.154 nm. The average distance between its parallel planes is 0.140 nm, which is much smaller than that 
in graphite. Single crystal diamond has a dense 3D4C structure comprising four nearest-neighbor polyacetylene- 
like zigzag chains or two nearest-neighbor twisted graphite layers 2 . Under static pressure, the highly crystalline 
varieties of graphite can be transformed to diamond via slipping, buckling, and cross-linking of the carbon 
sheets 3 " 9 with an sp 2 — > sp 3 bonding transition. During the last two or three decades, tremendous attention has 
been focused on synthesizing a variety of new carbon allotropes, including 0D3C fullerenes 10 , 1D3C nanotubes 11 , 
and 2D3C graphene 12 in all-sp 2 networks. Novel fullerene structures have been discovered in natural coal 13 , rock 14 
and detonation soot 15 . 

In addition to these well characterized carbon allotropes, there exist a number of proposed and synthesized 
carbon structures that require further exploration. Recently, we identified two 3D3C chiral framework structures 
in all-sp 2 bonding networks comprising 3-fold and 4-fold helical chains connected by ethene-type (H 2 C=CH 2 ) 
^-conjugation 16 . These structures contain matching helical chains of complementary chirality with ethene-type 
planar n- conjugation, and they are energetically more favorable than previously proposed structures with helical 
chains of same chirality (e.g., K 4 17 ) or zigzag chains (e.g., H-6 18 and bct-4 19 ) with twisted n bonds. This finding 
suggests a new approach to constructing covalent network structures. In particular, the matching helical chain 
configuration offers a new strategy for solving structures of yet unidentified carbon phases seen in recent 
experiments, including carbon blacks, soot, or similar materials 20 " 30 . 

Here we report by ab initio total-energy and phonon calculations 31 " 35 a new type of 3D3C crystalline carbon in 
an all-sp 2 network that comprises 6-fold helical carbon chains with alternating single and double bonds, which are 
distinct from the polyyne-like alternating single and triple carbon-carbon bonds in linear sp-carbyne 36 . This 
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Figure 1 | Chiral crystalline modification of carbon in an all-sp 2 3D3C bonding network in R3m {D\ d ) symmetry. Top (a) and side (b) view of rh6 
carbon, which comprises three 6- fold helices or three zigzag benzene rings as its building blocks. In hexagonal representation its lattice parameters are a = 
6.902 A, c = 3.470 A, and carbon atoms occupy the 18 h (0.8805, 0.1195, 0.5576) position. Symbols R and 6 represent right-handed helices and closed 
hexagon, respectively. It topologically corresponds to a (3 X 3 X 1) superlattice of graphite in AA stacking as shown in (c). 



structure can be regarded as twisted graphite in AA stacking consist- 
ing of zigzag hexagonal benzene rings connected by an ethene-type 
planar n- conjugation like as in the 3-fold and 4-fold carbene 16 . 
Phonon and electronic band calculations show that this new carbon 
phase is dynamically stable and is a semiconductor with a band gap of 
0.47 eV, in contrast to the semimetallic nature of graphite. Simulated 
x-ray diffraction (XRD) patterns show an excellent match between 
this newly identified carbon structure and the condensed carbon 
phase discovered in detonation soot 28 " 30 . These results establish a 
new carbon phase and offer insights into its outstanding structural 
and electronic properties. 

Results 

We first present structural characterization of the new carbon phase, 
which has a 6-atom rhombohedral primitive unit cell with equilib- 
rium lattice parameters a = 4.1491 A, a = 112.55°, occupying the 
6 h (0.4381, 0.7965, 0.4381) position, thus termed rh6 carbon (or 6- 
fold carbene). In hexagonal representation, it has an 18 atom hexa- 
gonal unit cell (see Fig. la and lb) with lattice parameters a = 
6.902 A, c = 3.470 A, occupying the 18 h (0.8805, 0.1195, 0.5576) 
position. It has a polymerized carbyne-like structure containing three 
6-fold helical carbon chains; however, instead of having the polyyne- 
like alternating sp-type single and triple carbon-carbon bonds as 
those in linear carbyne 36 , the right-handed (or left-handed) helical 
chains in rh6 carbon have an alternating sp 2 -type single and double 
carbon-carbon bonds (see Fig. la). The sp — » sp 2 bonding state trans- 
ition produces a large energy gain of about 0.46 eV per carbon atom 
(see Fig. 2a). Meanwhile, three right-handed (or left-handed) helical 
carbon chains bond together forming a zigzag benzene ring. 
Therefore, rh6 carbon also can be regarded as a twisted graphite in 



AA stacking consisting of three zigzag hexagonal carbon rings with a 
(3X3X1) superlattice of graphite (see Fig. lc). However, in contrast 
to the uniform bond length of 1.42 A in graphite, there are two 
distinct carbon-carbon bond lengths, a longer bond of 1 .483 A assoc- 
iated with a C(sp 2 )-C(sp 2 ) single bond in the zigzag benzene rings 
and a shorter bond of 1.359 A associated with an ethene-type 
C(sp 2 )=C(sp 2 ) double bond 16 between the zigzag benzene rings 
(these bonds are different from that in 3 -fold and 4-fold carbene 
with single bonds along the chains and double bonds between the 
carbon chains 16 ). Also, there are two different bond angles, 113.45° 
for ZC-C-C in the zigzag benzene rings, and 123.16° for ZC-C=C out 
of the zigzag benzene rings. These structural parameters on average 
are similar to those in graphite, but they are closer to those in 3 -fold 
and 4-fold carbene 16 , thus rh6 carbon can also be considered a 6-fold 
chiral carbene. 

We have calculated the total energy of rh6 carbon as a function of 
volume (see Fig. 2a) in comparison with the results for diamond and 
several other carbon structures that are all in an all-sp 2 bonding 
environment, including K 4 17 , H-6 18 , bct-4 19 , cR6 16 , cT8 16 , and graph- 
ite. The energies for fee C 60 fuller ene 37 and (5,5) carbon nanotube are 
also shown for comparison. Graphite has all its 2p z orbitals perfectly 
aligned, leading to a strong aromatic n bonding interaction between 
neighboring carbon atoms, which makes it the most stable carbon 
crystal. Conversely, all the 2p z orbitals are misaligned in the K 4 car- 
bon, making it the most unfavorable crystalline carbon form 16 . 
Meanwhile, rh6 carbon has one-third of its bonds with well-aligned 
2p z orbitals between the zigzag benzene rings and, as a result, it is less 
stable than diamond, graphite, C 60 fullerene, and (5,5) nanotube, but 
is almost as stable as the 4-fold cT8 carbene and more favorable than 
3 -fold cR6 carbene with the same R3m symmetry. 




Figure 2 | Energetics and kinetics of selected carbon structures, (a) Energy versus volume for iC 4 , H-6, bct-4, cR6, cT8 and rh6 in 3D3C bonding networks 
compared to those of 1D2C carbyne, 2D3C graphite and 3D4C diamond. The energies for fee C 60 fullerene and (5,5) carbon nanotube are also shown for 
comparison, (b) Enthalpy per atom (relative to that of graphite) for rh6, rh6-II, and diamond versus pressure. The insert picture is rh6-II structure. 
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Table I | Calculated equilibrium structural parameters (space group, volume Vq, lattice parameters a and c, bond lengths 

d c -c), total 

energy E fof , bulk modulus B 0 , and electronic band gap E g for cR6 carbon, cT8 carbon, rh6 carbon, graphite, and diamond at zero 
pressure, compared to available experimental data 38 



Structure 


Method 


V 0 (A 3 /atom) 


a (A) c (A) 


d c -c 


(A) 


Etot (eV) 


Bo (GPa) 


Eg (eV) 


Diamond (Fd3m) 


GGA 


5.60 


3.552 


1.538 




-9.018 


451 


5.36 


Exp 38 


5.67 


3.567 


1.544 






446 


5.47 


cR6 (R3m) 


GGA 


8.78 


7.122 3.597 


1.352, 


1.490 


-8.502 


268 


2.95 


cT8 (14 Cornell 


GGA 


8.40 


5.937 3.808 


1.351, 


1.488 


-8.551 


283 


2.41 


Rh6 (R3m) 


GGA 


7.96 


6.902 3.470 


1.359, 


1.483 


-8.550 


299 


0.47 


Rh6-ll (R3m) 


GGA 


5.94 


7.012 2.509 


1 .452, 


1.536 


-8.454 


408 


1.04 


Graphite (P6s/mmc) 


GGA 


8.81 


2.462 6.710 


1.422 




-9.045 


280 




Exp 38 


8.78 


2.460 6.704 


1.420 






286 





It should be noted that there is a dense form of rh6 (named rh6-II 
as shown in Fig. 2) with lattice parameters a = 7.012 A, c = 2.509 A, 
occupying the 18 h (0.7495, 0.1252, 0.4582) position. The rh6-II 
phase can be obtained from rh6 phase via local bond rotation (Fig. 
SI in Supplementary Information) and stabilized under pressure 
above 12 GPa (see Fig. 2b). It is even more stable than graphite above 
70 GPa due to the sp 3 -like bonding nature in rh6-II state. However, 
this structural transformation rh6 — » rh6-II is reversible (Fig. S2 in 
Supplementary Information), and upon decompression rh6 phase is 
recovered. The calculated equilibrium structural parameters, total 
energy, bulk modulus, and electronic band gap for rh6, rh6-II, cR6 
and cT8 are listed in Table I, compared to available experimental 
data 38 . 

To understand the mechanical stability of rh6 carbon, the corres- 
ponding elastic constants Cn, C 33 , C 44 , Ci 2 , and C i3 are estimated to 
be 630, 90, 128, 288, and 139 GPa, respectively. It clearly meets the 
mechanical stability criteria given by C 44 > 0, C n > |C 12 |, and 
(C u + 2Ci 2 )C 33 >2Cj 3 for hexagonal phase 39 . We further examined 
dynamic stability of rh6 carbon by calculating the phonon dispersion 
curves. No imaginary frequencies were observed throughout the 
entire Brillioun zone (results at 0 GPa are shown in Fig. 3a), con- 
firming the dynamic stability of rh6 carbon. The highest phonon 
frequency 1611 cm" 1 for rh6 carbon is almost equal to the highest 
phonon frequency of 1610 cm" 1 for graphite 40 , which reflects the 
graphite-like structural bonding character of rh6 carbon. 

The increased bonding connectivity in the 3D3C rh6 carbon is 
expected to open an electronic band gap in contrast to the semime- 
tallic state in the 2D3C graphite structure. A hybrid density func- 
tional method based on the Heyd-Scuseria-Ernzerhof scheme 



(HSE06) 34 has been used to calculate electronic properties. It is seen 
that the calculated HSE06 band gap (5.36 eV) for diamond is very 
close to the experimental data 38 of 5.47 eV, suggesting the validity of 
HSE06 method in predicting band gaps for diamond. Our calculated 
electronic band structure of rh6 carbon (see Fig. 3b) indeed shows a 
band gap of 0.47 eV with the conduction band bottom and valence 
band top located along the T- A and A-H direction, respectively. This 
moderate band gap makes rh6 carbon highly desirable for a wide 
range of applications in electronic, optical, and energy conversion 
devices. 

We now establish the experimental connection of rh6 carbon. 
Several experiments in recent years have discovered a new carbon 
structure in detonation soot or similar materials 20 " 30 , but have not 
been able to identify its structural features that hold the key to further 
understanding and exploration of this new carbon phase. 
Experimental x-ray diffraction (XRD) of a TNT/dissel oil detonation 
soot 28 (see Fig. 4b) revealed the presence of a considerable amount of 
amorphous carbon and several crystalline phases: the strongest peak 
around 26° comes from the graphite (002) diffraction; a weak peak 
around 43.6° is attributed to the diamond (111) diffraction; a prom- 
inent sharp diffraction peak at 30°, however, cannot be assigned to 
any known carbon phase such as graphite, diamond or fullerenes 25 " 27 . 
A similar sharp diffraction peak at 30° was also found in the chimney 
soot 29 and the trinitrotoluene/cyclomethylenetrinitramine detona- 
tion nanoparticles 30 . The high intensity and sharpness of this unex- 
plained XRD peak suggests that a new carbon phase has been 
consistently produced in these experiments. In Fig. 4a, we show 
simulated XRD patterns for graphite, diamond, rh6 carbon, cT8, 
cR6, H-6, bct-4, K 4 and fee C 60 described in Fig. 2a, and compare 




rAH KTMLH r A H KTM LH 



Figure 3 | Phonon and electronic band structures, (a) Calculated phonon band structure of rh6 carbon at zero GPa. (b) Calculated electronic band 
structure of rh6 carbon at zero GPa. 
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Figure 4 | X-ray diffraction (XRD) patterns, (a) Simulated XRD patterns 
for graphite, diamond, rh6 carbon, and the other structures cT8, cR6, H-6, 
bct-4, K 4y and fee C 60 . (b) Experimental XRD patterns for TNT/diesel oil 
detonation soot 28 . X-ray wavelength is 1.5406 A with a copper source. 

with the experimental XRD pattern 28 (see Fig. 4b). Our simulated 
XRD patterns for graphite and diamond are consistent with the 
corresponding experimental XRD peaks. Most importantly, the main 
peak for rh6 carbon matches almost perfectly the XRD peak at 30° for 
the newly identified carbon phase. Beside of rh6 carbon, it is shown 
that the main peak of 30.8° for K 4 carbon is also close to 30°; however, 
the K 4 carbon is dynamically unstable 41 . These results suggest that the 
new rh6 carbon is among the likely candidates for the unidentified 
carbon phases found in recent detonation experiments 28 . 

It is well-known that Raman spectroscopy is also an appropriate 
method to study the characteristics of new phases of carbon. To 
provide more information and characters for possible experimental 
observation, we also simulated the Raman spectra of rh6 carbon and 
compared the results with different sp 2 carbon structures (see Fig. S3 
in Supplementary Information). We find that the Raman spectrum 
of rh6 carbon presents a main peak A lg at 1605 cm" 1 and a weaker 
shoulder peak E g at 1580 cm" 1 instead of the E 2g mode at 1585 cm" 1 
in graphite 42 . These features may be helpful for identifying the new 
carbon phases in experiments 43 . 

Discussion 

In summary, we have identified a new type of 3D3C crystalline 
carbon in an a\\-sp 2 network in R3m (D\ d ) symmetry. It consists of 
three 6-fold carbon helices connected by three zigzag hexagonal 
carbon rings in hexagonal representation, and can be regarded as a 
chiral crystalline modification of graphite. Phonon and electronic 
band calculations indicate that this new carbon phase is dynamically 
stable and is a semiconductor with a band gap of 0.47 eV. Simulated 
XRD patterns provide an excellent match between this rh6 carbon 



phase and the previously unexplained carbon structure discovered in 
recent detonation experiments 28 " 30 . These results suggest that the 
possible presence of rh6 phase in soot carbon as well as the ultrafme 
diamond, graphite, and amorphous 28 . 

Recently, a one-step, gas-phase, catalyst-free detonation of hydro- 
carbon (C 2 H 2 ) method was developed to produce gram quantities of 
pristine graphene nanosheets 44 . It suggests an effective approach to 
synthesizing rh6 carbon via thermal decomposition or detonation of 
carbonaceous materials such as aromatic and polymeric hydrocar- 
bon molecules. Successful production of rh6 carbon in large quantity 
is expected to lead to greatly increased research interest and activity, a 
situation reminiscent of the early days of the discovery and ensuing 
study of fuller enes. 

Methods 

Calculations were carried out using the density functional theory as implemented in 
the Vienna ab initio simulation package (VASP) 31 . The generalized gradient 
approximation (GGA) developed by Armiento-Mattsson (AM05) 32 were adopted for 
the exchange-correlation potential. The all-electron projector augmented wave 
(PAW) method 33 was adopted with 2s 2 2p 2 treated as valence electrons. A plane-wave 
basis set with a large energy cutoff of 800 eV was used. Forces on the ions are 
calculated through the Hellmann-Feynman theorem allowing a full geometry 
optimization. The energy minimization is done over the atomic and electronic 
degrees of freedom using the conjugate gradient iterative technique. Convergence 
criteria employed for both the electronic self- consistent relaxation and the ionic 
relaxation were set to 10~ 8 eV and 0.01 eV/A for energy and force, respectively. A 
hybrid density functional method based on the Heyd-Scuseria-Ernzerhof scheme 
(HSE06) 34 has been used to calculate electronic properties. Phonon dispersion curves 
and Raman spectrum are calculated using the package MedeA 35 with the forces 
calculated from VASP. The phase conversion barrier was calculated using a gener- 
alized solid-state nudged elastic band method 45,46 with the cell and atomic position 
optimization. 
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